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Abstract 

We study the classical generalized g[„ Landau-Lifshitz (L-L) model with special 
boundary conditions that preserve integrability. We explicitly derive the first non-trivial 
local integral of motion, which corresponds to the boundary Hamiltonian for the SI2 L-L 
model. Novel expressions of the modified Lax pairs associated to the integrals of motion 
are also extracted. The relevant equations of motion with the corresponding boundary 
conditions are determined. Dynamical integrable boundary conditions are also examined 
within this spirit. Then the generalized isotropic and anisotropic Landau-Lifshitz 
models are considered, and novel expressions of the boundary Hamiltonians and the 
relevant equations of motion and boundary conditions are derived. 
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1 Introduction 



Numerous investigations have been devoted to the issue of incorporating non-trivial con- 
ditions that preserve integrabihty both in discrete [1]-[10], and continuum integrable 
systems [1], [11]-[18]. The central purpose of the present article is the study of classical 
integrable models when general boundaries that preserve integrabihty are implemented. 
Among the various classes of integrable models we choose to consider here a prototype 
model, that is the generalized classical continuum Heisenberg or Landau-Lifshitz (L-L) 
model. This model may be thought of as the immediate classical analogue of the XXX 
(XYZ for the anisotropic case) quantum spin chain [19, 20], whereas higher rank gen- 
eralizations may be seen as continuum limits of known high rank quantum spin chains. 
Although much attention has been devoted to the investigation of the quantum models 
with integrable boundary conditions not much progress has been made -from the al- 
gebraic point of view- on their classical continuum counterparts. Here, we consider the 
classical continuum case, we identify the boundary Hamiltonian, and the relevant bound- 
ary Lax pair, for the SI2 L-L model, utilizing primarily the algebra that rules the model, 
that is the classical reflection algebra. Our study concerns not only typical c-number 
reflection matrices, but also dynamical reflection matrices, which give rise to dynamical 
type boundary conditions for the aforementioned model. The generalized gl^ L-L models 
are also examined within this spirit. 

It is worth noting that the significance of the particular study stems primarily from 
the fact that it provides novel results for a wide class of classical integrable models 
associated to the g[„ algebra. Note that such anisotropic (trigonometric) models may 
be appropriately mapped to A^^li affine Toda field theories (see e.g. [21]). Moreover, 
the present investigation provides a first systematic description on the issue of integrable 
continuum limits of discrete integrable models, that contain boundary type terms or 
other distinct local terms, such as the ones arising also in the case of integrable defects. 

The outline of the article is as follows: in the next section we briefly review the model 
with periodic boundary conditions, as well as the relevant fundamental ingredients (see 
also [21]). The Lax pair formulation, and the construction of the associated integrals of 
motion through the associated algebras are reviewed. In section 3 we review Sklyanin's 
generic algebraic frame [1] describing classical models with boundaries that preserve 
integrabihty. Based on this framework we explicitly derive the associated Hamiltonian 
with suitable integrable boundary terms recovering also some of the expressions presented 
in [1]. Note that in [1] only diagonal boundary terms were treated, whereas here the most 
general boundary terms that preserve integrabihty are derived. Dynamical boundaries 
are also examined within this context. Note that our results are consistent with the 
classical continuum limits of the relevant quantum discrete Hamiltonians. In section 4 
we review the construction of modified Lax pairs in the presence of integrable boundaries 
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discussed in [2]. Relying on this framework, and using generic solutions of the reflection 
equation -c-number and dynamical- we are able to derive the Lax pair associated to 
the extracted Hamiltonians. This way the consistency of the whole procedure is fully 
ensured. In section 5 the continuum limit of the XXZ open spin chain is considered 
leading to the boundary antitropic sl2 L-L model. This further ensures the validity of 
the continuum limit process followed. We then examine the isotropic and anisotropic qI^ 
L-L models as continuum limits of the 0[„ and llg(gl„) open spin chains respectively, and 
obtain the classical continuum Hamiltonians, and the associated equations of motion and 
boundary conditions. 



2 The isotropic 5(2 Landau-Lifshitz model 

Let us briefly review the continuous isotropic Landau-Lifshitz (L-L) model with peri- 
odic boundary conditions, associated to the 5(2 classical algebra (see also [21]). After 
introducing the basic ingredients of the model and setting up our notations, we recall 
the Lax pair formulation for the classical integrable Hamiltonian system, and discuss the 
systematic means of constructing the whole tower of local integrals in involution. 

The physical quantities of the model are described by vector-valued functions S{x) = 
(5*1 (x), 5*2 (x), 5*3 (x)) taking values on the unit 2-sphere 

S'{x) = j2sKx) = l- (2.1) 

Note that throughout the text we shall also use the following combinations of Si{x) 

S^ix) = ^iS,ix)±tS2ix)). (2.2) 

The equations of motion associated to the isotropic Landau-Lifshitz model, which is our 
main interest here, are of the form: 

The fields Si{x) obey boundary conditions, which are taken to be such that either Si{x) 
become periodic, i.e. Si{x + 2L) = Si{x), or consider the fields and their derivatives to 
be zero at the endpoints (Schwartz boundary conditions). 

The Poisson structure of the phase space for the physical quantities Si{x) is given by 
the Poisson brackets 

{Sa{x), S,{y)} = 2ieabcSc{x)S{x - y), (2.4) 
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where Sabc is the totally antisymmetric Levi-Civita tensor with value £123 = 1. The 
Hamiltonian of the model is given by 

1 f I fdSiY fdS2V fdSs^ 



^^-jy i^iarj n^rj nari 1^^- p-^* 

The equations of motion in the Hamiltonian form are expressed as 

^ = {H, S). (2.6) 

Other physical integrals of motion include the momentum, which is given by 

P = / ' 7 dx, (2.7) 



1 + ^3 

and the total spin of the model in the case where periodic boundaries are considered. 



2.1 The Lax pair formulation 

Within the Lax pair formulation of a classical integrable Hamiltonian system one first 
defines the auxiliary linear differential problem, which reads as 

= U(x,t,A)^(x,t) 
^M/(a;,t) = Y{x,t,X)^{x,t). (2.8) 

In general, U and V are n x n matrices. Their entries contain dynamical fields, their 
derivatives and possibly, the spectral parameter A. The compatibility condition of these 
two equations leads to the so-called zero curvature condition 

dt\]-d,Y+[\], V] = 0, (2.9) 

which provides the equations of motion of the system under consideration. 
One then constructs the monodromy matrix 



T{x, y, A) = Fexp (^j^ l]{z)dzj , (2.10) 

being a solution of the equation (2.8). Assume that U obeys the classical linear Poisson 
algebraic relation [21] 

{U,(x,A), Ufe(2/,/i)} = Kfe(A-/i), Ua(x,A)+Ufc(y,^)]5(x-2/), (2.11) 

then it follows that the monodromy matrix satisfies the quadratic algebraic relation, 

{Ta{x, y, t. A), Tb(x, y, t, /i)} = [vabiX - /i), Ta{x, y, t, X)Tb{x, y, t, /i)]. (2.12) 
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rab is the classical r-matrix corresponding to the Hamiltonian system, and satisfies the 
classical Yang-Baxter equation [22]. The conserved charges may be obtained via the 
expansion of t(A) = trT{X) in powers of the spectral parameter, A. It can also be shown 
via (2.12) that these charges are in involution, that is they satisfy [21] 

{t(A), t(/i)} = 0. (2.13) 
In the case of the L-L model, the classical r-matrix has the simple form [23] 

r(A) = ^, (2.14) 



where V is the permutation operator: V{d ® &) = & ® a. 

We shall restrict ourselves for the moment in the case of the SI2 L-L model. In this 
case, the permutation and Lax operators are respectively 
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(2.15) 



We also note here the explicit form of the V-operator 

^^"^ = 2P'^-2Aa^'^- 



(2.16) 



Inserting the Lax pair operators (2.15) and (2.16) into the zero curvature condition (2.9) 
yields exactly the equations of motion (2.3). 



2.2 Integrals of motion 

As already mentioned, the local integrals of motion may be extracted through the ex- 
pansion of the trace of the monodromy matrix in powers of the spectral parameter A. A 
key property of the monodromy matrix that is crucial in what follows is that is satisfies: 

^T(x, y, A) = U(x, A)T(x, y, A). (2.17) 

Let us consider the following ansatz for the monodromy matrix: 

T{x, y, A) = (1 + W{x, A))e^(^'^'^)(l + W{y, X))-\ (2.18) 

W and Z are purely off-diagonal and diagonal matrices respectively. We also assume 
that W, Z are expressed as: 

00 00 
W{x,X) = J2^''Wn{x), Z{x,y,X)= J2 A"^n(a:,y). (2.19) 

n=0 n=—l 
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Our main aim henceforth is to identify the elements Wn, Zn, and hence the integrals 
of motion. It is technically convenient to split the Lax operator into a diagonal and an 
off-diagonal part as 

Substituting the ansatz (2.18) into the relation (2.17), and splitting the resulting equation 
into a diagonal and an off-diagonal part one obtains 

dW 

— + WVd - UdW + WUaW -Ua = 0, 
ax 

^=Ud + lJaW. (2.21) 

Plugging in the explicit expressions for Ua,Ud, we end up with the following relation for 
W 

^ + tWSs (^^ - TiS'a^ + S+a-) + ^WiS-a^ + S-^a-)W = 0, (2.22) 
ax X X X 

with a^, (T^ being the familiar 2x2 Pauli matrices: 

^-C:)' "^(oi)' -^C:)' '^-^^^ 

Inserting the expansion (2.19) in the equation above, one determines all W^s. We are 
only interested in the first three terms of the expansion, which are given by: 



Oil/X) : Wo 
C(A°) : -- 



-a\ 1-^3 2S+ 



a J 2S- I + S3 

-a'' 
-a' 



^ ^ ' a'-ia'^S" / U ' ^ ^ 



The first three terms suffice in order to compute the first two integrals of motion, namely 
the momentum and the Hamiltonian. To complete the computation we should also 
determine Z, provided by the following equation 

dZ 

— = IJ, + 1],W. (2.25) 

Substituting the relevant expansion of W and the explicit forms of U^, Ua into the equa- 
tion above we conclude 



Z{x, y,X) = —ix- y)a, + V A""^ / {S+a^ + 



a+)Wn{z)dz. (2.26) 
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By using the exact expressions for Wn, one may finally determine Z, order by order. 
In particular, as noted above, at order C(l/A) we found 

Z^, = ^{x-y)a\ (2.27) 

Here we set x = L and y = —L. The explicit form of Z_i is important since it indicates 
the leading contribution of as A — > 0, a result which will be greatly used below and in 
the forthcoming sections. Moving on to the next orders, one naturally arrives at the first 
two integrals of motion, as expected. More specifically, by working with the conventions 
defined above one concludes that 

(Zo)ii oc P, and {Zi)u oc H. (2.28) 

Extra care is needed while computing the integrals of motion above, taking into account 
the Schwartz boundary conditions at the endpoints for the fields and their derivatives, 
namely {Si{±L), Sl{±L)) — 0. Note that in the case of non-trivial boundaries extra 
terms, emerging from the bulk part, should be taken into account, given that the fields 
do not vanish in that case at the endpoints. This will be transparent in the subsequent 
sections. 



3 Implementing integrable boundaries 

We shall now briefly describe the relevant machinery needed to implement special bound- 
aries into a classical model in a way that integrability is ensured. After reviewing the 
general setting we focus on the particular example of interest, that is the Landau-Lifshitz 
model with integrable boundaries. 

3.1 Algebraic setting 

We review Sklyanin's formulation [1] in order to exhibit the integrability for classical 
models with special boundaries. Suppose that we have the monodromy matrix of the 
model at hand, endowed with the Poisson structure (2.12). Let also K'^i^X) be c-number 
(non-dynamical) representations of the classical reflection algebra [1, 11], 

0= [ri2(Ai-A2), i^(Ai) i^2(A2)] 

+i^i(Ai) ri2(Ai + A2) i^2(A2) - i^2(A2) ri2(Ai + A2) i^i(Ai). (3.1) 

Our notation is such that K+{X) = K{-X, k") and i^"(A) = K{X, k+); k± 
are generic free boundary parameters (see also below in the text). One may then define 
a modified transition matrix as [1] 

r{x,y,X)=T{x,y,X) K-{X) f{x,y,X), (3.2) 
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where T(A) = T A). The modified monodromy matrix satisfies the classical version 
of the refiection equation [1, 11]: 



{ri(Ai), r2(A2)} = [ri2(Ai - A2), r(Ai) r2(A2)] 

+ri(Ai) ri2(Ai + A2) r2(A2) - r2(A2) ri2(Ai + A2) ri(Ai). 



(3.3) 



The generalized transfer matrix reads as 



t{x,y,\) = tr{K+{\) Tix,y,\)}, 



(3.4) 



and it is shown to satisfy [1] 



{t{x,y,Xi), t{x,y,\2)} = 0, 



Ai,A2 G C, 



(3.5) 



hence it may be interpreted as the generating functional of the conserved integrals of 
motion. 

By adopting the ansatz (2.18) and setting x = 0, y = —L for the boundary points, 
the generating functional of the local integrals of motion takes the form 



In tr{(l + iy(0))-^i^+(A)(l + W^(0))e^(°'"^)(l + iy(-L))-ii^-(A)(l + H^(-L))e-^(°'-^)}, 



where W are the same as before, but with A — )■ —A. It is this expression that one expands 
in powers of A, in order to derive the local integrals of motion. 

3.2 The boundary Hamiltonian 

We now proceed in deriving the integrals of motion for the Landau-Lifshitz model with 
integrable boundaries. We shall be using the following c-number representation of the 
classical reflection algebra [3] 



As in the periodic case to derive the integrals of motion, one expands the generic object 
(3.6) in powers of A. We shall only present here the flnal results, whereas the technical 
details of the derivation can be found in the Appendix A. One should keep in mind that 
extra terms, emerging directly from the bulk, cancel out suitably some purely boundary 
terms providing eventually a quite simple boundary contribution to the Hamiltonian. 

In the case of open boundary conditions, the flrst integral of motion becomes trivial. 
This is actually expected, given that this is essentially the momentum, which is not a 



In tr 



{i^+(A) T(0,-L,A) K-{\) f (0,-L,A)} 



(3.6) 




(3.7) 
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conserved quantity anymore. The second integral of motion, i.e. the Hamiltonian is 
computed to be (see also [1], where only diagonal boundary terms are considered) 

^1 = -\ £^ (^^^ dx - ^(2k-5i(-L) - S,i-L)) + ^(2k+5i(0) - S3(0)).(3.8) 

The first term in the expression above is just the bulk term, while the rest are total 
boundary contributions. 



3.2.1 Dynamical boundary conditions 

We shall now discuss the case of dynamical degrees of freedom, attached at the ends of 
the system. To achieve this we shall consider the dynamical solution of the reflection 
equation [1]: 

K(A) = L(A) K{X) L-^(-A), (3.9) 



where we define 
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L = AI+|^g2_^ (3-10) 

The L matrix satisfies the quadratic algebra (2.12) with the Yangian r-matrix. The 
elements S^, S"*^ apparently satisfy the classical 5I2 algebra. Note that special limits of 
the generic L matrix lead to the Discrete-Self- Trapping (DST) model, or the Toda model 
(see e.g. [17] and references therein), so this way one treats a generic class of dynamical 
boundaries. 

We shall consider here for simplicity, but without loss of generality, the dynamical 
boundary attached to the left end of the system. The right end of the system will be 
described by the trivial reflection matrix oc I. Hence, the generalized transfer matrix 
will be of the form: 

t{X) = tr, [UX) K~{X) T-i(-A)] . (3.11) 
It will be convenient for the following computations to express the K-matrix as: 



K-(A) ocI + AB + C(A2), 1= , (3.12) 



where we define 




X = -4§3 + hsl - 1 - 4k"§3(§+ + 
Y = 8§+ + — f8k-(§+)2 - 2k-§2 _ 4§3§^ 

Z = 8S- + ^ (8k"(§")2-2k"§^-4§3§~). (3.13) 
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It is clear that the elements X, Y, Z contain the dynamical degrees of freedom attached 
to the boundary. Expanding appropriately the generalized transfer matrix as in the 
previous section, we end up to the following Hamiltonian with distinct dynamical terms 
attached to one boundary: 

I, = --ji — \ dx + -S^-L)Z+-S~i-L)Y--Ssi-L)X. (3.14) 

Note that such boundary terms would have emerged for the other end of the system as 
well, after implementing a similar right dynamical reflection matrix. More precisely, a 
dynamical matrix for the right boundary would lead to extra boundary terms in 
expressions (3.14) at x = of exactly the same from of as the ones at x = —L, but with 
r^-r, k-^k+ in (3.13). 



4 The modified Lax pair 

When integrable boundary conditions are implemented, the Lax pairs associated to the 
integrals of motion are accordingly modified. The systematic construction of the modified 
Lax pairs was presented in [2]. In what follows, we briefly review the results of [2], and 
then apply the formalism in the case of the s[2 Landau-Lifshitz model with integrable 
boundaries. 



4.1 Reviewing the construction 



Recall first the construction of the V-operator associated to a given integral of motion 
for a classical integrable model with periodic boundary conditions. Using (2.12) one 
formulates the following Poisson structure: 



Ta{L,-L,X), Ub(x,/i) 
where we define 



dM{x,X,fi) 
dx 



+ 



M{x,L,-L,X,fi), Ufe(x,/i) 



(4.1) 



M(x, A,/i) = Ta{L,x,\)rab{X - /i) Ta{x,-L,\). 



(4.2) 



More details on the derivation of the latter formula can be found in [21]. Recalling now 
that t{\) = trT{X) it naturally follows from (4.1) and (2.9) that 



V(x,A,/i), U(x,A) 



with 



V(x, A,/i) =t \\) tra(Ta{L,x,\) rab{X,fi) Ta{x,-L,\)y 



(4.3) 
(4.4) 
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In the case of open boundary conditions one may prove that a generahzed Poisson 
structure holds [2], i.e. 



{ra(0,-L,A), Ub(x,/i) 
where we now define 



^(x, A,/i)+ M<j(x,A,/i), Ub(x, /i) 



M(x,A,/i) = T{0,x,X)rab{\- f-i)T{x,-L,\)K~{\)f{0,-L,X) 
+ T(0, -L, \)K~iX)fix, -L, \)rabiX + /i)r(0, x, A). 



Finally, bearing in mind the definition of t(A), and (4.5) we conclude 

V(x, A, /i), U(x, /i) 



In t(A), U(x,/.)^ = ^^^^^ + 



where 



(9x 



V(x,A,/i) = t-i(A) trJi^+(A) M,(x,A,/i)). 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



This is the explicit form for the V-operator in the case of generic integrable boundary 
conditions. One expands V in powers of A in order to obtain the modified operator 
associated to each integral of motion of the model under consideration. 



4.2 Modified Lax pairs for the L-L model 

We are now in the position to determine the boundary Lax pair for boundary L-L model. 
The classical r-matrix associated to the L-L model is proportional to the permutation 
operator, and tVaVat = I, then V can be expressed in a simple form as 

V(a;,A,/i) = —^T{x,-L,X)K^{X)T~\0,-L,-X)K+{X)T{0,x,X) 
A — jj, 

+ ^—^T-\0,x,-X)K+{X)T{0,-L,X)K-{X)T-\x,-L,-X). (4.9) 
X + n 

The latter expression is valid for all classical models associated to the Yangian classical 
r-matrix, proportional to the permutation operator. 

Explicit computation shows that the V-operator for any point x ^ 0, —L reduces to 
the familiar bulk operator (2.16). In any case, we are mostly interested in computing the 
V-operator exactly at the boundary points, that is Xb = (0, —L). We shall only present 
the final results here, and postpone the heavy technical details of the computation until 
the Appendix C. At the end points the V-operator has the following form (we have 
multiplied the result of the expansion with |) 

Vfe(xb) =V(xb)+(5V(xb), (4.10) 



11 



where V is the bulk operator (2.16) and 
6Y{-L) -- 
6Y{0) -- 



1 


dx 


2 


2^ 




1 


dx 


2 


2j2 





\r{S+{-L) - S-{-L)) -\rS^{-L)-S- 
k-53(-L) + -\r{S+{-L) - 

k+(5+(0)-S-(0)) -k+S3(0) -5-(0) 
k+S3(0) + S+(0) -k+(5+(0)-S-(0)). 



-L) 



(4.11) 



From the zero curvature condition, and by requiring 5Y = (see more details on this 
argument in [2]), we obtain the equations of motion described in (2.3), and the non-trivial 
boundary conditions (see also [1] for only diagonal boundary conditions): 



Si 



dx 
dSi 

dx 
dS2 
dx 



Si 
Si 

S2 



dS2 

' dx 
dS, 

dx 
dSi 

dx 



-S2{-L) 
' -Sii-L) 



2k- 



■Ssi-L) 



2k' 



■S2{-L). 



(4.12) 



Of course one may easily check that the same equations of motion, and boundary con- 
ditions are extracted from the Hamiltonian (3.8) via (2.6). Note that in obtaining the 
boundary conditions we took into account that the derivative of the Casimir with respect 
to X is zero. Note that similar equations of motion are obtained for the other end of the 
system (a; = 0), but are omitted here for brevity. The entailed boundary conditions are 
as expected mixed ones. 



4.2.1 Dynamical boundary conditions 

The explicit computation of the modified V-operator in this case follows exactly the 
previous section's computations via the expression (4.9), so the result is quite straight- 
forward, as long as we keep in mind that the classical dynamical reflection matrix is 
now expressed as in (3.12). Recall that we restrict our attention here to one boundary 
Xb = —L, then the final expression for the boundary operator is given as 



Vfe(-L) = V(-L) +5V(-L), 



(4.13) 



where we define 



dx 



ZS+{-L) - YS~{-L) -ZSsi-L) - 2XS~{-L) 
YSsi-L) + 2X5+(-L) -ZS+{-L) + YS~{-L) 



(4.14) 



where X, Y, Z are defined in (3.13). In this case the relevant boundary conditions, 
entailed from the conditions 6Y = 0, read as 



dS, 
dx 



S- 



dS2 
dx 



i{Z - Y) 



Ss{-L)+XS2{-L) 



12 



S: 



Si 



dx dx 

oi— 02 



-XSi{-L) 



Y + Z 



S^i-L) 



dx 



dx 



Y + Z 



^2 



-L) - (4.5) 



Needless to mention that these boundary conditions emerge also from the dynamical 
Hamiltonian (3.14) through (2.6). 



5 Integrable continuum limit 

We shall describe here a systematic means of obtaining classical continuum limits of 
quantum discrete theories for generic boundary conditions along the lines discussed in 
[20]. Assume a collection of operators assembled in matrices Lu, acting on "quantum" 
Hilbert spaces labeled by i and encapsulated in a matrix "acting" on the auxiliary space 
Vi. For any quantum space q they obey the quadratic exchange algebra [24, 25] 

R12 Liq L2q = L2q Lig R12 , (5-1) 

where operators acting on different quantum spaces commute, and R satisfies the Yang- 
Baxter equation. The form of the monodromy matrix T is then deduced from the co- 
module structure of the YB algebra 

Ta = LaN La2 ■ ■ ■ Lai , (5-2) 

and thus naturally obeys the same quadratic exchange algebra (5.1). 
First consider that the R matrix has a classical limit as 

R=l + hr + 0{h'^), (5.3) 

with r satisfying the classical Yang-Baxter equation (see also e.g. [22, 20]). We may now 
establish that T has a classical limit by considering in addition the classical counterpart of 
L, which then satisfies the quadratic Poisson algebra emerging directly as a semi-classical 
limit of (5.1), after setting j:[A, B] — > {A, B}. It reads 

{La{Xi), L,(A2)} = KbiXi - A2), LaiXi) L,(A2)] . (5.4) 

The classical discrete monodromy matrix is apparently of the same form as in (5.2). The 
exchange algebra for T takes the form 

{Ta,n} = [rab, TaTb] . (5.5) 

This quadratic Poisson structure implies that the traces of powers of the monodromy 
matrix tr(T'^) generate Poisson-commuting quantities identified as classically integrable 
Hamiltonians. 
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Now that we have discussed the classical limit we may proceed to the continuum limit 
of discrete theories with open boundary conditions. In this case the modified monodromy 
matrix has the form 

nX)=n\) K^{\)T~\-X), (5.6) 

where T is given by (5.2), T satisfies the classical refiection equation (3.3) and K~ is a c- 
number solution of the refiection equation (3.1). Introduce a suitable spacing parameter 
6 : 0{S) 0{j^). Let us also express the L matrix as 

Lan{\) = l + Sl]a„i\)+5^Ui^{X) + ... 

L~l{-\) = l-<5U..(-A) + 52t/£)(-A) + ... (5.7) 
It then naturally follows for the monodromy matrix and its inverse 

r(A) = l + 55^U,„(A) + 52 5^U,„(A)U,„(A) + 52^t/£)(A) + ... 

n n>m n 

T-\-\) = i-5^U„„(-A) + 52 5^U,„(-A)U,„(-A) + 52^f/i^,H-A) + ... 

n 

(5.8) 



For the following we consider that 



(5.9) 



Based on the latter formulas it is clear that terms with powers of 6 bigger than the number 
of summations go to zero in the continuum limit. This is the so-called "power counting" 
argument presented in more detail in [20]. So the continuum limit of the monodromy 
matrix (5.8) becomes 

T(0,-L,A) = Pexp{y" dxl]{x)}. (5.10) 

We conclude that the continuum limit of the discrete modified monodromy matrix reduces 
to the continuum analogue of (5.6), with T given in (5.10). The open transfer matrix, 
the generating function of the charges in involution as usual is 



t(A)=tr,(i^+(A)r,(A)), (5.11) 



where a c-number solution of the classical refiection equation. Having said these it 
is clear that continuum limits of discrete Hamiltonians would provide legitimate Hamil- 
tonians of integrable continuum theories. In the two examples below we make this clear 
comparing also with the results of the the two previous sections. 



14 



5.1 The open XXX chain 



It is quite straightforward now to check that the expression of the boundary Hamiltonian 
of the L-L model may be directly extracted from the quantum XXX open chain as an 
appropriated continuum limit (see e.g. [19, 20]). Recall the open XXX Hamiltonian [1], 
(see Appendix B for more details on the derivation): 



7V-1 



J2 [^"^^u+^y.+i+^y.^i 



2e 



2k-((T++(Ti-)-CT^ 



2^H 



Nj 



'N 



(5.12) 

The bulk term, after introducing appropriate coherent states see e.g. [19, 20], and after 
making the following identifications: 



and 



cr 



Six + S), 



(5.13) 



reduces to the bulk part of the Hamiltonian (3.8), whereas the boundary contributions of 
the left and right boundaries reduce exactly to the boundary terms of the Hamiltonian 
(3.8). Note that an implicit rescaling A — )■ 6~^X takes place. It is worth stressing here that 
continuum limits of discrete systems with boundary terms should be taken with particular 
care in order to have a sensible result. This is an intriguing issue also encountered in the 
case of integrable systems with local defects. 



5.1.1 Dynamical boundary conditions 

As in the non-dynamical case the classical Hamiltonian may be also directly derived as 
a continuum limit of the open XXX chain with a left dynamical boundary term (see also 
Appendix B): 

1 ^""^ 1111 
where we define 

X = -4(§3 + 1) + -^(2^1 + 1 - 2k-{§-, §3} _ 2k-{§+, §3}) 
Y = 8§+ + (8k-(§+)2 - 2k" (§2 - 1) _ 2{§3, 
Z = 8§" + ^(8k-(§-)2 - 2k-(§2 - 1) - 2{§3, 

D = —(^2k" (§++§") -§3), (5.15) 

where { , } denotes the usual anti-commutator. Notice that the quantum dynamical 
K-matrix is again of the form (3.9), however we considered for convenience the quantum 
Lax operator to be of the same structure as in (3.10), but with an additional | term. The 
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quadratic Casimir is chosen to be zero, so that the results are compatible with the classical 
continuum case. This means that we choose to consider the spin zero -non-compact- 
representation of s[2- Note also that we considered here a trivial right boundary oc I, 
a dynamical matrix for the right boundary would lead to extra boundary terms in 
expressions (5.14) at x = of exactly the same from of as the ones at a; = — L, but with 
-^+, k- ^ k+ in (5.15). 

The boundary term of the quantum Hamiltonian above reduces to the boundary term 
of the classical continuum Hamiltonian, ensuring the consistency of the whole process 
with integrable continuum limits. There are some extra terms, which can be seen as 
"quantum corrections" , and at the classical limit they vanish. 



6 Generalized boundary L-L models 
6.1 The boundary anisotropic 5(2 L-L model 

This section serves mostly as a further check on the consistency of the continuum limit 
process we considered in the previous section. In particular, we shall start with the open 
XXZ Hamiltonian with generic boundary terms, and then we shall appropriately take 
the continuum limit in order to obtain the corresponding continuum Hamiltonian. 

Let us first recall the associated i?-matrix: 

R(X) = ( ^'""^^^ + ^ + 2) sinh(^i)(T- A 
^ ' \ sinh(/i)a+ sinh(A - ^ + f )/ ' 

The generic reflection matrices, associated to left and right boundaries, [12], are given in 
Appendix B. 

The Hamiltonian is then defined in Appendix B (see eq. (B.6)), and in a more explicit 
form is expressed as: 

« = ^ E (<'!''?« + "'"j^^ + )<'^J«) + a^hj^-) ( - »sh(»r)<^J + 2k-af) 

Now recalling the identifications (5.13), the arguments presented in [20], and setting 

sinh(/i) = Sy/J + cosh(/i) = l-S^^ + ... (6.3) 

we may take the continuum limit of the latter expression and obtain the following classical 
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Hamiltonian (see also [1]): 



dx 



L 



dS_ 

dx 



JSl 



2 smh(i^ ) 
-cosh(ze+)53(0) + 2k+5i(0) 



- cosh{iC)S3{-L) + 2k"5i(-L) 



(6.4) 



2sinh(i^+) 

Utilizing the exchange relations among the classical spin variables we end up to the 
typical bulk equations of motion 



Si 

S3 = iiSi 



I S, 



dx'^ 
d^Si 

dx'^ 
d^S,. 



S. 



S, 



s, 



dx"^ 
d'S3 

dx'^ 
d^Si 



iJS2S3 

iJS^Si 



dx"^ dx^ 

and the associated mixed boundary conditions (see also [1]): 



(6.5) 



dSs ^ dS2 
'^2— J3 



S3 
Si 



dx 
dSi 
dx 
dS2 
dx 



-Si 
-S2 



dx 
dS3 
dx 
dSi 
dx 



cosh(i^ )a/J 



sinh(i^ 
cosh(z,^~)vO' 



S2{-L) 



-L sinh(z,^ 

'''^-S2i-L) 



Si{-L) 



2k- 



sinh(z^ 



-S3{-L) 



(6.6) 



=-L sinh(z,^- 

It is clear that in the isotropic limit we recover the equations of motion, and the associated 
boundary conditions derived in the previous section. Apparently analogous boundary 
conditions emerge for x = 0, which are omitted here for brevity. 

Having checked once more the consistency of the continuum limit process, we may 
now proceed with the 0[„ generalizations of the isotropic and anisotropic (trigonometric) 
L-L model. 



6.2 The boundary isotropic L-L model 

We may now generalize our analysis for the isotropic 0[„ model. The easiest way to 
obtain the desired results is to start with the Hamiltonian of the corresponding quantum 
spin chain, and take the suitable continuum limit. This is an absolutely legitimate 
way of extracting continuum integrable Hamiltonians as has been transparent from the 
description presented in the preceding sections. 

Let us first introduce the R and K matrices associated to the g[„ model. The 0[„ R 
matrix is given as 

n 

R(X) = X + V where ^ = ^ ^ij ® e^j, (6.7) 
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where we define: {eij)ki = SikSji. The K"^ matrices in this case are of the generic form 
(see also 



K~{X, r, k") =I + AB~, B" = ^(^-qen -c^e„„ + 2k-(ei„ + e„i) + cl 

K+{\, k+) = I - (A + 1)B+, B+ = ^ - c+en - c+e„„ + 2k+(ei„ + e„i) + cl 
where cf = + 1, ^ = - 1 c = 4k^^ + 1. (6.8) 

Then the boundary quantum discrete Hamiltonian is given as (see also Appendix B): 

n N-1 
i,j=l m=l 

- ^( - c^eSf - 4ei-) + 2k-(eS;) + ei^))). (6.9) 
Consider now the following identifications (see also [20]), 

eS^) ^ h,{x), eS7+') k,{x + 6) (6.10) 
where the elements kj satisfy the classical g[„ algebra: 

\^hj{x), Ikiiy)^ = (sjjkix) - Sjkhi{x)^5{x - y). (6.11) 
We then conclude about the continuum limit of the expression (6.9): 

1 " 1 / 

Xi = -- / rfx E^*i(^)^M^) + -cr/ii(-^)-C2"U(-i^) + 2k-(/i„(-L) + /„i(-L)) 

- ^ ( - c^^ii(O) - ctLniO) + 2k+(/i„(0) + /„i(0))) (6.12) 

where the prime denotes derivative with respect to x. 

From the latter expression the equations of motion and the corresponding generic 
mixed boundary conditions for the model are readily entailed: 

n 

4z(x) = E {l';i{x)l,,{x) - ll,{x)Ux)) (6.13) 
i=i 

and 



1 

E [U-L)l,,{-L) - l',^{-L)l,i{-L)) = - [ - cr (5n/u.(-L) - 5,,ln{-L)) 

(6.14) 

Analogous boundary conditions are clearly obtained for x = 0. 
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6.3 The boundary anisotropic L-L model 

We shall finally describe the anisotropic (trigonometric) L-L model starting from the 
open quantum spin chain. The i?-matrix associated to the spin chain under 
consideration is given by [26]: 

R{X) =p(^sinh(A + /i) I + sinh(A) u), (6.15) 

where V is defined in (6.7), U provides a representation of the Hecke algebra (see also 
[26]), and is defined as 

n 

\J= J2 (e^i ® e,, - g-^^"('-^\,, ® e,,) (6.16) 

where q = e^. We consider the reflection matrix written in terms of the representation 
of the "boundary" element e of the affine Hecke algebra (see e.g. [8]) 

K^{X) =x^{X)I + y^{X) e^, (6.17) 

where we define 

x~(A) = —6q cosh(2A + /i) — k~ cosh(2A) — 2 cosh(22C^) sinh(/i), 
y-{X) = 2sinh(2A) sinh(/i) 

x+(A) = -6^ cosh(-2A - /i) - cosh(-2A - 2^) - 2 cosh(2<+) sinh(/i), 
y+{X) = -2 sinh(2A + 2/i) sinh(/i) 

= -iQ^ + iQ^r'), = QiQ^r' + q^'Q^, (6.18) 

and Q^, are free boundary parameters. Also, choose a particular representation 
provided by 

= -{Q^y^en - Q^Cnn + ei„, + e„i. (6.19) 

The quantum spin chain Hamiltonian may be written then in terms of the affine Hecke 
elements as (see also [8]): 

Af-l 

H = ^ U, ,+1 + C- + C+ e+ , (6.20) 

i=i 

where = =F Q±-(Q±^-'i'-'^2^Lh(2tC='=) " ^S^^^ considering the identifications (6.10), and 
expressing U as [20] 

n 

U = P- I + ^(^1- g""3"(^-^'))e,, ® ejj (6.21) 

and 

/i = Sa, q~'9n{i-j) _ I _ gg^^^ _ ^^^^ ^ (g 22) 
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we obtain the following continuum boundary Hamiltonian, which may be seen as a de- 
formation of the isotropic g[„ model presented in the previous section: 



Xi = -- / dx(^ l[j{x)l'ji{x) + lii{x)ljj{x) + 2a ^ {ln{x)l'jj{x) - ljj{x%{x)) 

•^-^ i, j=i i^j=i i<j=i 

+ C-[- (g-)-^/ll(-L) - Q~lnn{-L) + lU-L) + lnl{-L)) 

+ - (g+)-^/ii(o) - g+u(o) + /i„(o) + /„,i(o)), (6.23) 



where = =F Q±_(Qjz)-f"2cosh(2ic^) " -^^^^ latter expression and the gt^ exchange 
relations the associated equations of motion are entailed 

n 

+ 2a XI - 2a X l'jj{x)lki{x) - 2a ^ l'^.{x)lik{x) + 2a ^ l'^^{x)lki{x), 

j<l j<k j>l j>k 

(6.24) 

with corresponding boundary conditions 

n 

J2 [U-L)l,u{~L) - l'^^{-L)l,i{-L)) = -aYhA-L)kk{-L) 

+a Y lni-L)lkii-L) + a ^ ljj{-L)lik{-L) - aj^ 
j<k j>i j>k 

-C 



-Q ^{5illik{-L) - 5iklii{-L)^ - Q^SniLki-L) - 5nklnl{-Lj^ 
SiiLki-L) - 5nkhi{-L) + 6nihki-L) - 5iklni{-L) , (6.25) 



and similarly for the other end of the theory at x = 0. With this we conclude our analysis 
on the generalized g[„ boundary L-L models. Note that similar generalizations can be 
applied in a straightforward manner in the elliptic case, but are omitted here for brevity. 

7 Boundary symmetries 

We focus here mainly on the isotropic 0[„ case, and briefly discuss the continuum ana- 
logues of earlier woks on boundary symmetries (see e.g [8]), although we have to mention 
that this is a whole separate subject of interest. 

We shall extract below the so-called boundary non-local charges which are realizations 
of the underlying classical reflections algebra. We shall first consider oc I, and K~ 
provided by the generic expression: 

K-{\) = \ + M. (7.1) 
A 
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The L and L ^ operators are expressed as 

Lo/(-A) = l + ^Po, + ^p2. + ... (7.2) 

We now consider the modified monodromy matrix T, and by expanding in powers of 
i we extract the boundary non-local charges (see also e.g. [8]). Let us start with the 
discrete T 

ToiX) = Lon{X) . . . Loi{X) K-{X) L^,\-X) . . . L-^(-A), (7.3) 
the expansion of the latter leads to 

r„(A) = rr + ir„<" + ir« + ^^. 

= Bo + ^(5^Po^Bo + lo$^Poi + 

i i 

+ p ( 5Z PoiPo,-Bo + Bo 5^ PoiPoi + Bo J] P^, + ^ PoiBoPo, + P 

i>j i<j i i,j j 

(7.4) 

To obtain the continuum limit of the latter expression we recall formulas (5.9) as well as 
the "power counting" argument, then 

/•O 

5^Po. ^ / rfxPo(x), 

(5' ^Po,Po, ^ I dxdy Po(a;)Po(y). (7.5) 
The modified monodromy matrix as well as the B, P matrices may be expressed as 

n 

T = GkiTki 

k,l=l 
n 

P = ^CfciPfci 

k,l=l 
n 

B = ^CfciBfe,. (7.6) 

k,l=l 

The entries Tki of the matrix are the non-local charges, which are realizations of the 
classical reflection algebra (3.3), (see also [27]). In particular, the continuum non-local 
charges emerging from (7.4) via (7.5) become: 

'kl - ^fci, 



n „0 pO 
'Tip = J2 dx¥km{x)M^l + J2^km dxFml{x), 

m=l-^~^ m=l ~^ 
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'kl 



^ / dx dy¥k^{x)¥^p{y)Mpi+ ^Mkm dx dy ¥rnp{x)¥.pi{y) 

n,p=l J^>V m,p=l "^^<2/ 

n „ 

+ J2 dx dyFkUx)Mmprpi{y), ... (7.7) 

Consider now the continuum expression for T with T given in (5.10), and take into 
account that for the L-L model 

U{x) = \f{x), (7.8) 
A 

where P(x) = J2ab=i ^ab^ab{x), and the entries fab{x) satisfy the classical gl„ algebra. It 
is then clear that expansion of the continuum T in powers of ^ would exactly lead to 
expressions (7.7) (see also the relevant discussion in section 5). Once more the consistency 
of the continuum limit is manifest. 

We may now show that Tj^P , which form a closed algebra themselves, provide an 
exact symmetry of the continuum open transfer matrix (see also [5, 8]), i.e. 

t(A), r,«] = 0. (7.9) 

We be gin our proof by considering the classical reflection algebra (3.3) as Ai — >■ oo: 

{B<, + i-ri^) + ..., 7^(A2)} = j^^l^ab, (Ba + ^rw + ...)r5(A2)" 

VabTbiXi) 

(7.10) 



Ai 


-A2 




1 


Ai 


+ A2 




1 


Ai 


+ A2 



Keep only the first order terms ^ in the latter formula (set also A2 = A), and bear in 
mind that [B^, 7^] = then, 

{%^^\ TbiX)} = Vab^aTb{\) - ^aTb{X)Vab + ^aVabTb{\) - Tb{\)Vab^a 

{r«, tn{rb{x))} = . . . =0 ^ {ri/\ t(A)} = o (7.11) 

and this concludes our proof on the exact symmetry of the transfer matrix. Note that 
the indices a, b refer to the auxiliary spaces, whereas the indices k, I denote entries of 
the matrix according to (7.6). 

In fact, the above algebra, formed by 7^'"/\ is the full 0[„ algebra for K~ oc I [8], and 
it breaks down to suitable subalgebras depending on the choice of K^. For instance, if 
we choose to be of the form: 

K{X) = dzag( a{X),...,a{X), 6(A), 6(A) V 
V"' ^ ' ^ V 

I n—l 

a(A) = -A + 2e, b{X) = X + i^ (7.12) 
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then the exact symmetry reduces to gli ® 0l„_/ [5, 8]. 

Finally, a more generic choice of matrix would further reduce the exact symmetry 
of the transfer matrix, these issues have nevertheless been examined in e.g. [5, 8], and we 
shall not further discuss them here. Note also that the symmetry for anisotropic models 
is discussed in detail in [8]. 

8 Discussion 

Let us summarize the main findings of the present investigation: we have been able 
to explicitly derive expressions of the first classical integrals of motion for a prototype 
integrable model with non-trivial boundaries, that is the 5(2 isotropic Landau- Lifshitz 
model. We considered both c- number (non-dynamical), as well as dynamical classical 
reflection matrices, which give rise to non-trivial, but still integrable, boundary terms in 
the Hamiltonian of the model. Compatibility of our results with the classical limit of 
the corresponding quantum discrete model Hamiltonian, i.e. the XXX open spin chain, 
further ensures the validity of our findings. In addition to the integrals of motion we have 
been able to derive novel expressions for the associated boundary Lax pairs for the SI2 
isotropic L-L model, following the prescription introduced in [2]. As expected both the 
boundary Hamiltonians as well as the associated Lax pairs lead to the same equations of 
motion and boundary conditions, verifying the consistency of the methodology followed. 

Having checked the consistency of the continuum limits of the XXX and XXZ open 
spin chains, leading to the isotropic and anisotropic SI2 L-L model respectively, we then 
generalized our analysis to the isotropic and anisotropic g[„ L-L models. More precisely, 
starting from the corresponding open spin chain Hamiltonians, with generic boundary 
terms, we considered the corresponding continuum limit, and extracted the associated 
classical boundary Hamiltonians as well as the equations of motion with the relevant 
boundary conditions. Finally, a brief discussion on the boundary symmetries is also 
presented. It is shown that the presence of special boundary terms suitably breaks 
the symmetry of the models under consideration, as also happens in discrete integrable 
models (see e.g. [5, 8, 27]). 
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A Derivation of local integrals of motion 



We provide here the main technical points on the derivation of the local integrals of 
motion. Recalling that the leading contribution of and as A — )■ comes from the 
Z\\ and Z\\ entries, we consider the following expansions 



;i + W^(0, A))-^J^+(A)(l + 1^(0, A))] ^^ = Y, A"/in, 

n=0 

oo 

;i + W{-L, \)Y^K~{X){\ + W{-L, A))] = 5^ A"/i, 



n=0 



Z(0, -L, A) - Z(0, -L, A) = -L + V(l - (-)")A" / S- {W,,^^)^^dz. 

-Ill A — • ./„r. 



n=l 

Hence, the expansion of the modified transfer matrix is given by 

In tr{/s:+(A)T(0, -L, A)i^-(A)f (0, -L, A)} 
1 

A^ 



(A.l) 



^ oo „o / CO oo \ 

-L + ^(1 - (-)")A" / SYWn+i)2idz + In 5^ 5^ hJimX''^'^ 

71=1 "^"^ \n=Om=0 / 



• (A.2) 



However, we also need to expand the logarithm in the equation above. As will be clear 
below, the first terms ho and ho are proportional to unit. This fact enables us to write 
the logarithm as 



oo oo 



oo oo 



In 5^ XI ^n^-A"""" = In 1 + X(/l„ + /l„)A" + ^rJlm\ 



n+rn 



(A.3) 



n=l 



vn=0 m=0 

where the hn, hn that appear above have been rescaled as^ 



n=l m=l 



hn 1— ■ 



7 hn 
hn T-: 

ho 



(A.4) 



and expand thus the logarithm into powers of A as 

/ oo oo \ oo 

In 5^X/in/^^A"+'" =$^/nA^ 



(A.5) 



v,n=0 m=0 



n=l 



where /„ provide essentially the boundary contributions to the integrals of motion for 
the left and right boundary respectively. It is interesting to observe that the boundary 
contribution decouples into terms associated with the left and right boundaries, that is 



no mixmg occurs 



/i = /ii + /ii, f2 = -hhl + h2 - Ihl + h2 



^This amounts to adding the term ln(— ) to the expansion of the modified transfer matrix. 
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/s = Ihf - hih2 + h3 + Ihl - hih2 + hs, 

(A.6) 

Next, we show explicitly how to compute the h^s, that is the expansion of the first 
equation in (A.l). 

We now compute the inverse of (1 + W). Let 

oo 

(l + l^)-i = 1 + ^A"F„. (A.7) 

n=0 

By demanding that the object above is the inverse, 

(1 + Fo + AFi + X^F2 + ---){l + Wo + XWi + \^W2 + ■■■) = 1, (A.8) 

we determine the quantities 

Fo = -1^0(1 + W^o)-' 

Fi = -(1 + ^0)1^1(1 + 1^0)-' 

F2 = -[il + Fo)W2 + FiWi]{l + Wo)-\ ■■■ (A.9) 
Recall also that W{X) = W{-X), then 

W2n+1 ^ -W2n+1, Tl = 0, 1, 2, ■ ■ ■ . (A.IO) 

By denoting {1 + W)-^ = 1 + '£n=o ^"^n, we find: 

Fo = -Wo{l + Wo)-' 

A = +(l + Fo)W^i(l + W^o)"' 

F2 = -[il + Fo)W2-F,Wi^{l + Wo)-\ (A.ll) 
Substituting the explicit forms of the relevant quantities we arrive at 
ho = 

h = 2k+ S,{0) - 53(0) + ^e ( 't|°g^(of^ + ^3(0)) 

= 2k+ sm - sm + ^-^(^5^(0) + smsm - sms2{o)), (a.u) 

and 

ho = iC 



h, = -2k- s,i-L) + s,i-L) - ^^ (^^^rrlBr^ ^ ^^^""^0 

-2k- S,{-L) + Ss{-L) - —-L^(^^S',{-L) + S^{-L)S',{-L) - S[{-L)S2{-L)) 

(A.13) 



These functions are not the rescaled ones defined in (A. 4). It is clear from the expansions 
of the left and right boundary contributions that the following symmetry holds 

hn = {-irhn: e^r, k+^k-. (A.14) 
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B The open spin chain Hamiltonian 



The open Hamiltonian emerges as a first derivative of the associated quantum transfer 
matrix of the model (see also e.g. [1]). Note that the quantum transfer matrix for the 
boundary discrete model is of the same structure as the classical one (3.4). Let us first 
consider the isotropic case The final expression of the open quantum Hamiltonian is given 
by: 

noc^^H, + -—^ + ^^^^^^0) (B.l) 



2 d\ 



where we define for the isotropic (g[„) case: 

= 'Pjk- (B-2) 

Assume that the ]K--matrix, solution of the quantum reflection equation, has the follow- 
ing generic form 

K" ocI + AB + C>(A^), (B.3) 
it is then clear that the boundary contribution is essentially the B-matrix. 
In the XXX case we have: 

c-number representation 

dynamical representation, (B.4) 

where the elements X, Y, Z, D for the dynamical case are defined in (5.15). Also the 
right boundary is defined as: 

^ ' l^-2k+(A + l) + ^ ' 

For the anisotropic case the generic expression one gets for the Hamiltonian is: 

^ , sinh/i dK^jX) , troiK^iO)Hm) 

iJ oc ^ H,,^, + 2sinh(ze-) 6X a=o + tvoKHO) ^^'^^ 





^-1 2k- 




^2k- 1 


-x + 


D Z 


Y 


X + D) 



where we define 

rr _ dRjj+i{X) 
J^jj+l - I 3j+^ 

In the XXZ case in particular we choose: 



A=0 



(B.7) 



K~{X) 



sinh(-A + iS,-) k" sinh(2A) 
k" sinh(2A) sinh(A + 



K+(X) = I sinh(A + + ze+) -k+ sinh(2A + 2^) 
^' ' -k+ sinh(2A + 2/i) sinh(-A-/i + 2e+ 
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C The modified V-operator 



To compute the modified V-operator, we begin with (4.9). The monodromy matrix 
satisfies T{x, x, A) = T~^{x, x, A) = I, so the results are greatly simplified at the boundary 
points. 

Substituting the ansatz (2.18) into (4.9) and taking into account the explicit form of 
the transfer matrix presented in Appendix A, (A.l), one derives the exact expressions 
for the V-operator. 

For the boundary point X;, = we have: 

V,j = Xr l-^AaBi, + -^CmA , (C.l) 

where we have defined 

A = 1 + W{0), B = {1 + W{0))-^K+{X), 

D = {l + WiO))-\ C = K+{X){1 + W{0)), (C.2) 

and Xr is the quantity appearing in the first line of (A.l). It is straightforward to 
compute the expression for V, since each one of the quantities above is known. 

In a similar fashion, for the boundary point Xb = —L: 

where we have now defined 

A = K-{X){1 + W{-L)), B = {1 + W{-L)y\ 

D = {l + W{-L))-^K-{\), C = {l + W{-L)), (C.4) 

and Xi is now the quantity appearing in the second line of (A.l). 
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